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Abstract. We prove existence and smoothness of the density of the solution to a non- 
hnear stochastic heat equation on L^iO) (evaluated at fixed points in time and space), 
where O is an open bounded domain in R** with smooth boundary. The equation is 
driven by an additive Wiener noise and the nonlinear drift term is the superposition 
operator associated to a real function which is assumed to be (maximal) monotone, 
continuously differentiable, and growing not faster than a polynomial. The proof uses 
tools of the Malliavin calculus combined with methods coming from the theory of 
maximal monotone operators. 



1. Introduction 

Let O C M"^ be an open bounded domain, and consider the stochastic semihnear 
parabohc evolution equation on L^(0) 

du{t) - Au{t)dt + f{u{t))dt = r]u(t)dt + BdW{t), u{0) = uq, (1.1) 

where A is the Dirichlet Laplacian on L?'{0), is a standard cyhndrical Wiener process 
on L^{0), / : M — )• M is continuous, increasing and of polynomial growth, is a positive 
number, B is a deterministic bounded linear operator on L'^{0), and uq is an L?'{0)- 
valued random variable (precise assumptions are given in the next section) . As explained 
in Remark 12.11 the above equation covers the case where f — rj is replaced hy f + g, 
where / is as above and g is a globally Lipschitz function. 

Under regularity assumptions on the initial datum and the coefficients, one can prove, 
using monotonicity methods, that (jl.ip admits a unique mild solution with continuous 
paths (see e.g. [1] and references therein). The purpose of this paper is to study the 
existence and regularity of the density (with respect to Lebesgue measure) of the random 
variable u{t, x), with (t, x) s]0, T] x O. In particular, we show that a sufficient condition 
for the existence of the density is that / is of class with polynomially bounded 
derivative. We do not claim that this condition is sharp (as a matter of fact, we conjecture 
that it is not, by far), but we do hope that our methods could be the starting point for 
further developments. The proof relies on techniques of the Malliavin calculus and 
on a priori estimates for solutions of approximating equations, obtained replacing the 
nonlinear drift term by its Yosida regularization. Furthermore, we show that if / is of 
class C"* with polynomially bounded derivatives, the density becomes smoother, as it 
is natural to expect. The proof of this still relies on a priori estimates, and requires a 
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further regularization of the drift, as the Yosida approximation does not necessarily have 
bounded derivatives. 

There is a large literature on problems of existence and regularity of densities for 
solutions to parabolic SPDEs with Lipschitz non-linearities by means of the Malliavin 
calculus (see [TU [TBI [ISl HSl EI] and references therein). On the other hand, much 
less attention has been dedicated to SPDEs with non-Lipschitz coefficients: the first 
article (relying on Malliavin calculus) we are aware of is the article [20], where a nonlin- 
ear one-dimensional stochastic heat equation on [0, 1] with polynomially growing drift 
and diffusion coefficients is considered. Using techniques of the Malliavin calculus, the 
authors prove, under quite general conditions, including a very weak non-degeneracy 
assumption on the diffusion coefficient, the absolute continuity of the law of the solution 
evaluated at fixed points in time and space. We note that existence and uniqueness 
of solution for the stochastic heat equation on [0, 1] driven by an additive noise with a 
measurable drift having polynomial growth has been proved in [ID]. In [8j, the absolute 
continuity for the law of the solution to a one-dimensional stochastic heat equation with 
Holder continuous diffusion coefficient and linearly growing drift is proved with com- 
pletely different methods. In particular, the proofs involve the Euler approximation and 
techniques of harmonic analysis. 

As far as applications of Malliavin calculus to non-linear parabolic SPDEs is con- 
cerned, we should mention that the existence of the density for the stochastic Cahn- 
Hilliard equation with non-Lipschitz coefficients has been investigated in [3|. Further- 
more, the existence and smoothness of the density for the solution to the stochastic 
Burgers equation with globally Lipschitz coefficients have been considered in |12l I23j . 

The paper is organized as follows. In Section [21 we provide the definition of mild solu- 
tion to our equation (jl.ip . state a well-posedness result for this equation and summarize 
some properties of a regularized version of in terms of the Yosida approximations. 
Section[3]is devoted to establish the random field counterpart (i.e. d la Walsh [22]) of our 
equation. In Section [H we collect some auxiliary results needed for the main theorem's 
proof. These correspond to a version of the chain rule for Malliavin derivatives, some 
properties of time-dependent evolution operators and estimates for some regularizations 
of the drift coefficient. Eventually, in Section [SJ we state and prove the main result of 
the paper. Theorem 15. II For this, we first deal with the Malliavin differentiability of the 
solution to our equation and, secondly, we study the invertibility of the corresponding 
Malliavin matrix. 

Notation. We shall write a < 6 to mean that there exists a constant N such that 
a < Nb. To emphasize that the constant N depends on the parameters pi, ■ ■ ■ ,Pm, we 
shall also write a <pi,...,p^ b. We shall write sup to denote both the supremum and 
the essential supremum. For a Banach space E, we shall denote by L^(£^) the space of 
S-valued random variables ^ such that E||^||^ < oo. 

2. WeLL-POSEDNESS and APPROXIMATION 

Let O be an open bounded domain in W^. spaces over the domain O will be 
denoted without explicitly mentioning the domain, e.g. stands for L'^{0). The norm 
in LP will be denoted by || • and the inner product of by (•, •)2, unless otherwise 
stated. 
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On a given stochastic basis (f^, J-", F = {J~t)o<t<T,^), with T a fixed positive number, 
let us be given the following semilinear SPDE on L^: 

du{t) - Au{t) dt + /(u(t)) dt = r]u{t) dt + B dW{t), n(0) = uq, (2.1) 

where A is the Dirichlet Laplacian on L^, / : M — t- M is continuous, increasing, and such 
that |/(x)| < 1 + for some p > (we denote the evaluation operator associated to the 
function / by the same symbol), ry is a positive number, is a cylindrical Wiener process 
on generating the filtration F, i? : — )• is a linear and bounded operator, and uq 
is an J"o-™6asurable L^-valued random variable such that E||uo||2 < expressions 
involving random quantities are meant to hold P-a.s. if not otherwise specified. 

Remark 2.1. Instead of ()2.ip one may equivalently consider an SPDE of the type 

du{t) - Au{t) dt + f{u{t)) dt = B dW(t), u{0) = uq, 

where / is quasi-monotone, i.e. such that / + 77 is monotone for some r/ > 0. In fact, 
we may write / = (/ + ??) — 77 = / — ??, with / monotone, thus obtaining (|2.1|) with / 
replacing /. Furthermore, note that equations of the type 

du(t) - Au{t) dt + f{u{t)) dt + g{u{t)) dt = B dW(t), n(0) = no, 

with / monotone (or quasi-monotone) and g (globally) Lipschitz are just particular cases 
of (j2.ip . In fact, one has 

if{x) + g{x) - f{y) - g{y)) {x - y) > {g{x) - g{y)) {x-y)> -Hs-Hc-o,! I^; - y\^, 

i.e. f + g is quasi-monotone, since f + g + i] is monotone, with rj = \\g\\(jo,i- (Here || • ||^o,i 
stands for the Lipschitz norm). 

We shall work with the so-called mild solution, whose definition we recall. In the 
following we shall denote by A the realization of the Dirichlet Laplacian on different 
function spaces. The same convention we are going to use for the semigroup S{t) := e*^, 
t > 0, generated by A. 

Definition 2.2. An L^-valued adapted process u is a mild solution to equation (|2.1|) if 
f{u) E L^{[0,T] — >■ L^) and it satisfies the integral equation 

u{t) + [ S{t- s) {f{u{s)) - r]u{s)) ds = S{t)uo + [ S{t- s)B dW{s) (2.2) 
Jo Jo 

for all t € [0,r]. 

As is well-known, in order for the stochastic integral in ()2.2p to be a well-defined 
(Gaussian) L^-valued random variable it is necessary (and sufficient) to assume that 

/ Tr {S{s)BB*S{sy) ds < 00 VtG[0,r]. (2.3) 
Jo 

Let us also recall that this condition is weaker than requiring Q := BB* to be trace- 
class (cf. e.g. [SI Ch. 2]). We shall actually work under the following stronger standing 
assumption, which guarantees, as we are going to recall, that (|2.ip admits a unique mild 
solution with paths in a space of continuous functions. 
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Hypothesis 2.3. It holds 



E sup 

t<T 



S{t- s)BdW{s) 







_ < oo Vq > 1. 

C{0) 



Conditions on B and O implying that this hypothesis is fulfilled are extensively dis- 
cussed in the literature (cf. e.g. [H Ch. 6] and [5l Ch. 2]). 

In order to state the well-posedness and approximation results we need, let us fix some 
further notation: we denote by Cg, 1 < g < oo, the space of adapted processes u with 
values in C{0) such that 



n||c, := (Esuplln(t)II^^^J < oo. 



The following global well-posedness result holds true (see e.g. [6l Thm. 7.13 and Rmk. 
11.23], as well as [U Prop. 6.2.2], for a proof). 

Theorem 2.4. Let q>l and assume that uq G L'^(C(0)). Then equation (|2.ip admits 
a unique mild solution u (z Cq. 

Let us now introduce the regularized SPDE 

dux{t) - Aux{t)dt + fx{uxit))dt = r]ux{t)dt + BdWit), u{0) = uq, (2.4) 

where fx, A > 0, stands for the Yosida approximation of /, i.e. 

fx{x) ■=l{x- M^)), ^a(x) := (/ + Xfy\x). 

Recall that fx is increasing, Lipschitz continuous with Lipschitz constant bounded by 
1/A, and fx^f pointwise as A ^ 0. Moreover, Jx is a contraction and fx = f°J\ 
(see e.g. [2] for a detailed discussion of the Yosida approximation). It is then clear that 
Theorem 12 .41 applies also to (|2.4p . yielding the existence and uniqueness of a mild solution 
Ux € Cq. Furthermore, the following convergence result holds true (see [H Prop. 6.2.5]). 

Proposition 2.5. Let q > 1 and assume that uq S h^{C{0)). Then there exists a 
constant N , independent of A, such that 



Moreover, one has 



Esnp\\uxml^^^<N. 



limEsup||nA(t)-u(t)||^-^ =0, 



where u & Cq denotes the (unique) mild solution to equation ()2.ip . In particular, one 
has, for any {t, x) G [0, T] x O, 

lim E\ux{t,x) - u{t,x)\'^ = 0. 

A->0 
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3. An ALTERNATIVE EXPRESSION FOR THE SPDE (|2.ip 

The pathwise continuity in t and x of the solution to equation (j2.2|) guaranteed by 
Theorem 12.41 allows us to pass to the random field formulation of (j2.2|) . i.e. as in [25], in 
a sense made precise in Proposition 13.11 below . 

It is classical that S{t) is a kernel operator for all t > 0, i.e. there exists a function 
]0,T] X O X O 3 {t,x,y) ^ Gtix,y), with G L^{0 x O) for all t G]0,T], such 

that, for any < t < T, 



Jo 

and 

||5'(t)||l^oo = \\Gt{-,-)\\L°°{OxO)^ 

where || ■ ||i-i^oo stands for the — )■ L°° operator norm. Since the semigroup S is 
contracting in L°° i.e. ||<S'(t)/||Loo < ||/||loo for all t > 0, and it holds 

Gt{x,0)= [ Gtix,y)dy = [S{t)lo]{x), 

one has 



o 



supGt(x,0) = ||5(t)lo||Loo < 1. (3.1) 
xeo 

We have the following result, where we use the terminology introduced in [22] for sto- 
chastic integrals. 

Proposition 3.1. For q > 1, let uq G h'^{G{0)) and denote the unique mild solution in 
Cq to equation ()2.2|) by u. Setting u{t,x) := [u{t)]{x), one has, for any {t,x) g]0,T]x O, 

u{t,x)= Gt{x,y)uQ{y)dy + I / Gt^s{x,y){riu{s,y) - f{u{s,y))) dy ds 
JO JO JO 

+ [ [ Gt^s(.x,y)W{ds,dy), 
Jo JO 

where W stands for a martingale measure with covariance operator Q = BB* . Moreover, 

E ( sup \u{t,x)\'^ ) < cx), (3.2) 

\{t,x)€OT J 

where Or := [0,T] x O. 

Proof. It is readily seen that we only have to prove that, for any t G]0,r], 

f S{t-s)BdW{s)= f I Gt^s{-,y)W{ds,dy) (3.3) 
Jo Jo JO 

as L^-valued random variables. Let us assume (without loss of generality) that, formally, 

W{t) = Y,e^Wk{t), (3.4) 

fcGN 
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where {e^j^gN is a basis of and {wfclfcgN is a family of independent standard one- 
dimensional Brownian motions. Then we have, by the integral representation of S{-), 

[ S{t-s)BdW{s) = y2 [ S{t- s)Be''dwkis) 

= Y.I I Gt^s{:y)[Be%y)dydwu{s), 

where the series of ordinary Ito integrals are no longer formal by virtue of ()2.3p . 

As far as the stochastic integral on the right-hand side of ()3.3p is concerned, we 
notice (see e.g. [7]) that it may be understood as a stochastic integral with respect to a 
cylindrical Q- Wiener process {Wh{t), he L^,t>0} in the sense of e.g. [TJ]. Namely, 
by definition, the latter is a centered Gaussian family of random variables such that, for 
any h G L^, the process {Wh{t), t > 0} is a Brownian motion with variance t{Qh,h) 
and, for all s,t > and h,g £ L^, 

E{Whis)Wg{t)) = {sM){Qh,g)2. (3.5) 

By an innocuous abuse of notation, this cylindrical Q- Wiener process will be also denoted 
by W . Let us introduce the Hilbert space Lg, which we define as the completion of 
with respect to {h,g)j^2 := {Qh,g)2- Then one can define the (real-valued) stochastic 

integral with respect to W of any Lg-valued square integrable process as follows: let 
{e'^jkeN be a basis of Lq and X G L^($7 x [0,r] — > Lq) (see e.g. [3 Sec. 2]), and set 

r [ Xis,y)W{ds,dy) := ^ l\x{s),e^) L^^dW,,{s). 

Jo Jo fcGN"^0 

Moreover, the L^-valued cylindrical Wiener process W of (|3.4p determines a cylindrical 
Q- Wiener process W, with Q = BB*, as follows: for any t >0 and h G L^, set 

ft 



^h{t) ■■=y^ I {Be\h)2dwk{s). 



k& 

It is immediate that Whit) is a centered Gaussian random variable and 

e(^^ / {Be^,h)2dwk{s)y = t^{e^,B*h)l = t\\B*h\\l = t{BB*h,h)2 = t{Qh,h)2. 

k •'^ k 
In a completely similar way one verifies the covariance condition ()3.5p . 

In order to prove (|3.3p . one just needs to take := Q~^/'^e^ , where denotes 
the pseudo- inverse of Q^^^, whence one easily verifies that W^k{s) = Wk{s) and 



Gt^s{;y)[Be'']iy)dy= / Gt-s{;y)[Qe']{y) dy. 
lO JO 

We have thus proved the identities 

f S{t-s)BdW{s) = y^ f f Gt-s{;y)[Be'']{y)dwi,{s)= ff Gt-s{;y)W{ds,dy). 
Jo ^ Jo Jo Jo Jo 

Finally, the estimate ()3.2p is an immediate consequence of Theorem 12.41 □ 



DENSITIES FOR PARABOLIC SEMILINEAR SPDES 



7 



4. Auxiliary results 

We collect here some tools that we shall need in the next section. In particular, 
we give a version of the chain rule for the Malliavin derivative, where a (Malliavin) 
differentiable random variable is composed with an increasing function of polynomial 
growth. The result, also without the monotonicity assumption, is certainly well-known 
to experts (cf. e.g. [131 P-36]), but we include a proof for completeness and because 
most standard references contain only a proof for functions of class (cf. e.g. 
Prop. 1.2.3]). Moreover, we prove some estimates for evolution operators generated 
by time-dependent bounded perturbations of the Laplacian and for regularizations via 
Yosida approximations as well as via mollification. 

As usual, a function : R — )■ R is said to be of polynomial growth if there exists 
p G N such that |0(x)| <l + \x\P for ah x G R. We shall denote by Cpoi(R) the space 

of functions (j) G C'^(R) such that cf), (j)' , . . . , (f)^^^ are of polynomial growth. It is not 
difficult to see that, by the fundamental theorem of calculus, one can equivalently say 
that Cj„i(R) is the space of functions <j) € C'^(R) such that 0^*^^ is of polynomial growth, 
and also that it is the space of functions (j) G C'^(R) for which there exists p G N such 
that 

i0(x)i + i0'(x)i + ... + i#)(x)i <i + |xr 

4.1. Malliavin calculus. We shall repeatedly use the following strong- weak closability 
property of the Malliavin derivative (cf. e.g. [181 P- 78]). We use standard notation and 
terminology (see e.g. [T7]). 

Lemma 4.1. Let A; G N, p g]1,oo[. Assume that lim.„^oo = X in L^i^L) and 
sup„gN ll^nljo'^'P < Then X G J])^'^ and X„ — )> X weakly in H^'P . 

We shall denote the Gaussian Hilbert space "supporting" the Malliavin calculus by 
H. Moreover, we shah say that X G if X G B'^'P for aU p > 1. 

Lemma 4.2. Let X G 1])^'°°, / G Cj^oi(R) and increasing. Then f{X) G B^'^ and 
Df{X) = f{X)DX. 

Proof. Assume, without loss of generality, that |/(x)| -|- |/'(x)| ^ 1 + |x|p, with p G N. 
Let /a, a > 0, denote the Yosida approximation of /. Since fx^f pointwise and / is 
continuous, one has fx{X) f{X) P-a.s. as A ^ 0. Moreover, recalling that |/a| < |/|, 
hence that, for any g > 1, 

i/A(x)-/(x)r<2'?i/(x)r <i + ixr, 

and E|X|'^P < oo, the dominated convergence theorem yields 

limE|/A(X)-/(X) 1-^ = Vg>l. (4.1) 

Appealing to the inverse function theorem, it is easy to show that / G implies G Cb, 
and 



1 + A/'(Ja (!•))' ' l + A/'(Ji(i)) 
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In particular, since fx is Lipschitz continuous, one has fx £ C^, so that the "classi- 
cal" chain rule (see e.g. [HI Prop. 1.2.3]) implies that fx{X) G ©1'°° and Dfx{X) = 
f'^{X)DX . Let us check now that, for any (7 > 1, it holds 

supE||D/a(X)||^ <oo. (4.2) 

A>0 

By the above expression for and /' > it immediately follows that 

\f'x{x)\<\f'{Jx{x))\<l + \Jx{x)\P<l + \x\^, 

hence, by Cauchy-Schwarz' inequality, 

E||D/a(X)||^ < {^f'x{X)\^'if'^{^\DX\\^^Y''^< (1 + E|X|2P9)^/^(E||Z)X||^^)^/^ 

The right-hand side is finite and independent of A, thus (j4.2p is proved. In particular, 
dUD and imply that f{X) € B^'^ and Dfx{X) converges to Df{X) weakly in 

U[H) as A —7- 0, for all g > 1. Since the weak limit is unique, in order to prove that 
Df{X) = f'{X)DX, it suffices to show that 

limE\\Dfx{X) - f'{X)DX\\j, = 0. 

A — 

Observe that 

E\\Dfx{X) - f'{X)DX\\j, = E(|/;(X) - f'{X)\''\\DX\\j,) 

< {E\fiiX)-f'iX)f'>f\E\\DXf^)'^\ 

and that, by the expression of f'^, one has /j^ — t- /', hence fx{X) f'{X) P-a.s. as 
A — >■ 0. Moreover, taking into account that 

\f'x{x) - nx)\ = I ^ fxnMx)) " ^'^^^1 - + - 1 + 1^1"' 

one also has 

|/;(X)-/'(X)P'?<1 + |X|2P^ 

hence, by the dominated convergence theorem, f'x{X) f'{X) in L^"^ as A — t- 0, and the 
proof is finished. □ 

4.2. Time-dependent evolution operators. Using the notation introduced in Sec- 
tion [21 let F : [0, T] L^, and consider the following linear evolution equation on 

dy{t)-Ay{t) + F{t)y{t) = 0, y{s) = yo€L'', 0<s<t<T. (4.3) 

Here L^, p G [l,oo], denotes the set of (equivalence classes of) functions (j) G such 
that (j) > a.e.. The evolution operator U{t,s) is then defined by y{t) =: U{t,s)yQ. 

Proposition 4.3. For any < s < t < T , the following properties hold true: 

(i) U{t,s) is positivity preserving, i.e. yo > implies U{t,s)yQ > 0; 

(ii) U{t, s) < S{t — s), i.e. yo > implies U{t, s)yo < S{t — s)yo; 

(iii) U{t,s) is ultracontractive, i.e. its — > L°° norm is finite. 
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(iv) U{t,s) is a kernel operator, i.e. there exists a function k : [0, T]^ x 
such that 



[U{t,s)(p]{x) = k{t,s;x,y)(f>{y)dy. 



o 

Proof, (i) Let yo > and y{-) be a strong solution (without loss of generality) of ()4.3p . 
where we assume, for simplicity, s = 0. Taking the scalar product with y~ and integrat- 
ing with respect to time, we obtain, denoting the norm by || • ||, 

h\y~it)f+ [ \\Vy-{r)fdr+ [ {F{r)y-{r),y-{r)) dr = \\y^f, 
^ Jo Jo 

hence \\y~{t)\\^ < 2 Wy^W^ = 0, i.e. y{t) > for all t € [0,T]. (ii) Let z be the solution 
to z' — Az = 0, z(0) = yo- Then one has 

j^{y - z){t) - A{y - z){t) + F{t)y{t) = 0, {y - z)(0) = 0, 
hence, taking the scalar product with (y — z)"*" and integrating, 

\\{y{t) - z{t))+f + /* ||V(y(r) - z{r))+f dr + /* (F(r)y(r), (y(r) - z(r)) + > ds = 0, 
Jo Jo 

which yields, recalling that, by (i), y(r) > for all r, y(t) < z{t) for all t. (iii) Note that, 
by (i), one has 

\U{t,s)yo\ = \U{t,s)y+ - U{t,s)yQ\ < \U{t,s)y+\ + \U{t,s)yQ\ 

= U{t,s)y+ + U{t,s)yQ =U{t,s)\yo\, 

therefore 

\U{t,s)yo\<Uit,s)\yo\<Sit-s)\yo\GL^, 

which immediately implies U{t,s)yo € L°°. (iv) is a direct consequence of (iii), thanks 
to a classical criterion of Dunford and Pettis (see e.g. [Ill Ch. XI, §1]). □ 

4.3. Regularizations. In the next Lemma we use the notation introduced immediately 
after Theorem 12.41 



Lemma 4.4. Let f G C™[(]R) be an increasing function such that \f^'^\x)\ < 1 + \x\^ 
for all 77, = 0, 1, . . . , m. Then the following properties hold: 

(i) fx and Jx belong to C"^(M); 

(ii) for A < 1, there exists g € N, independent of X, such that |/{"^(a;)| < 1 + \x\'^ for 
all n = 0,1, ... ,m. 

(iii) /|"^ converges pointwise to /("^ as A — )■ 0. 

Proof. Since Jx = (I + A/)"^, and /' > 0, the inverse function theorem implies that, if / 
is of class C™, then also Jx is of class C". Moreover, the identity A/a = I — Jx implies 
that also fx is of class C™, hence (i) is proved. 

(ii) The polynomial growth of fx is obvious by the inequality |/a| < |/|. Moreover, 
recalling that fx = f°J\^ that Jx is of class C™ and is a contraction, we have 

1/1(^)1 = \f'{Jx{x))\ \j'x{x)\ < \njx{x))\ < 1 + \Jx{x)r < 1 + 
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Taking into account that Xf" = J" and that 

fx=f"{Jx){Jx? + f'iJx)Jx, 

one gets (1 + Xf{Jx))f'{ = f"{Jx){J'x)\ which imphes 

1/^(^)1 < I/"(^a(x))| < 1 + \Jx{xW < 1 + l^r. (4.4) 

Unfortunately it does not seem possible to extend such elementary arguments to obtain 
the polynomial growth of We can nonetheless argue as follows: by Arbogast's 

formulsQ (see e.g. [9]) one has 

where the sum is taken over all possible combinations of 62> • • • i &n} C N U {0} such 
that 

bi + 262 + • • • + nbn = n and bi + b2 + ■ • • + bn = k. 
In particular, note that there is only one possible term of the series containing f'{Jk), 
precisely the one corresponding to 61 = 62 = • • • = ^n-i = 0, 6„ = 1, that is f'{Jx)J^\ 
Similarly, there is only one possible term containing f^"'\Jx), precisely the one corre- 
sponding to 61 = n, 62 = &3 = • • • = &n = 0, that is /'^")(Ja)( Jj^)". Recalling that 
A/{"^ = we have 

(1 + A/'(Ja))/1") = /(")(Ja)(4)" + (4.5) 

hence 

|/1"^|<|/(")(Ja)| + |s„-i|, 

where Sn-i is a finite sum of terms involving only f^^\Jx) and powers of jf \ for 
k = 1, . . . , 71 — 1. Using once again that Xf^^ = —J^\ we conclude that Sn-i is a finite 
sum of terms involving only /^'^^(Ja) and powers of Xf^\ for k = 1, . . . ,n — 1. Taking 
A < 1, since |/('')(Ja(x))| < 1 + \x\p for all /e = 1, . . . , n, recalling that f" satisfies (lilil) . 

/ON 

we obtain that there exists 93 S N such that |/^ ^ 1 + \x\'^''^. By iteration one ends 

up with \f^\x)\ ^ 1 + \x\'^'' for each k = 0,1, . . . ,m. Since m is finite, this implies the 
claim, (iii) It is known that fx ^ f pointwise as A — )■ 0. That f'^ f pointwise as 
A ^ has been proved in Lemma Passing to the limit (pointwise) as A ^ in (j4.5p 
we obtain 

lim/f) = /(") + lim5„_i, 

A^O ^ A->-0 

where we have used that J\ ^ I and 

lim J'\{x) = lim . , , , = 1. 

A^o ' A^o 1 + Xf'{Jx{x)) 



This formula is better known as Faa di Bruno's formula. The latter attribution, however, seems to 
be historically incorrect. 



DENSITIES FOR PARABOLIC SEMILINEAR SPDES 11 

The claim is proved if we show that —?■ pointwise as A — )• 0. For n = 2 one has 

5*1 = 0, hence the claim holds. For n > 3, we observe that each term in ^^-i is of the 
form 

where c is a positive number, 1 < i,hi < n — 1, and /12, /13, . . . , ^n-i ai'e nonnegative 
integers with at least one of them, say hs, greater than 1. Let 2 < a := hs- Recall that 

J^^^ = -A/f \ 

therefore the generic term of 5^-1, hence Sn-i itself, converges to zero as A — )• 0. □ 

Let us introduce mollifiers in the following (standard) way: for C, € C°°(M) positive, 
with support contained in [—1, 1] and Jj^C = 1, set, for any /3 > 0, Cp{x) '■= /^~^C{^/ P)- 

Lemma 4.5. Let f : R ^ R be such that \f{x)\ < 1 + \x\p, and //? = /* C/3; /?< 1- 
Then one has 

\f^{x)\ < N{l + \x\P) VxGM, 
where the constant N does not depend on (5. 

Proof. Assume that |/(x)| < A^i(l + \x\^)- By the triangle inequality and the properties 
of Q one has 

\fp{x)\< [ \f{x-y)\C(s{y)dy<N, [ {i + \x - y\P)C(s{y) dy 
<iVi(l + |xn + iVi / \y\PC^{y)dy 
= N,{l + \xn + N,l3P [ \y\PC{y)dy 

jR 

<Ni{2 + \x\P). □ 

5. Existence and smoothness of the density 

In this section, we provide sufficient conditions on the data of the problem implying 
that the law of u{t, x) is absolutely continuous with respect to the Lebesgue measure 
and that the corresponding density is a differentiable function. More precisely, we will 
prove the following result. 

Theorem 5.1. Assume that uq € C{0) and let u € HgeN '^'J mild solution to 

(|2.1|) . Furthermore, assume that there exists 7 € (0,2) such that, for all x G O, there 
exists a constant > such that, for any t € (0, 1), 

g{x,t):= [ \\Gs{x,-)\\l, dsyc^t-y. (5.1) 

JO 

(a) If f £ CpQ[(IR) is increasing, then for any {t,x) €]0,T] x O, the random variable 
u{t, x) is absolutely continuous with respect to the Lebesgue measure. 

(b) Moreover, if for some integer m > 2, f G C™[(M), then the density of u{t,x) 
belongs to C^'^iR). 
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As we shall see below, the assumption (|5.ip is not needed to prove Malliavin regularity 
of u{t,x). It is instead needed to prove finiteness of negative moments of the Malliavin 
matrix (which reduces to a real random variable in the present setting). Moreover, as 
explained in Remark 1 5. 9 1 in order to prove the existence of density, condition (j5.1|) can be 
slightly weakened. Nevertheless, for the sake of conciseness and clarity, we have decided 
to state only one condition of the term g{x,t). 

Remark 5.2. The hypothesis on the initial datum in the previous theorem may be relaxed 
to uo € hi{C(0)) for ah ^ > 1, and uq G L°°{0 B^'°°). However, it does not seem 
natural to assume the initial datum to have such regularity. 

Let us give some examples of domains O and covariance operators Q = BB* satisfying 
condition (|5.ip above. 

Example 5.3. Let d = 1, O = (0, 1) and B = Id. Then ()5.ip holds with 7 = ^ (see e.g. 

m (A.3)]). 

Example 5.4. Let O = (0, vr)''. Define, for any k = {ki, . . . ,kd) e N'^, 

efc(x) := (^^^ sin(A;iXi) • • • sin{kdXd), x € O. 

Then, it is readily checked that the family {efcl^jgp^d is an orthonormal basis of L^(0) 
such that 

-Ack = \k\'^ ek, 

where \k\^ — kf + '-' + kj. Set B = {I - A)"™ for m > 0. Then, since Q = (I - A)-^™, 
one has 

g{x,t) = / / Gsix,y)[QGsix,-)]{y)dyds 
Jo Jo 

= I" Y.{l + \kfr^{Gs{x,-),ek)l2ds 



2 

kern 



Using the fact that |efc(x)| is uniformly bounded with respect to k and x, one easily 
verifies that m > ^ — 1 implies that the latter series is finite. Moreover, we have that 



.2,|fc|2 ^ 2t\k\^ ^ 2t\k\' 



l + 2t\k\^ - l + 2r|A;|2' 
Hence 

g{x,t) >t J] (1 + |A;|2)-™ (1 + 2T\k\Y' \ek{x)\ 



keNd 
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and this series can be bounded from below by any of its summands, such as the corre- 
sponding to /c = (1, . . . , 1) € N'^. Therefore, 

g{x,t)>c^t, with := (l + (i)~™(l + 2r(i)"^ 0y sin(xi)---sin(2;rf). 

Since x S (Oj-tt)"^, it is clear that Cx > and this implies that condition (|5.1|) is fulfilled 
with 7 = 1. 

Before turning to the study of the Malliavin differentiability of the solution to our 
equation, let us recall that the underlying Gaussian space on which to perform Malli- 
avin calculus is given by the isonormal Gaussian process on the Hilbert space H := 
L^([0,T] -^q)> which can be naturally associated to the cylindrical Wiener process 
W with covariance Q. With a slight (but harmless) abuse of notation we shall write W 
instead of W for notational convenience. 



5.1. Malliavin difTerent lability. The purpose of this subsection is to prove regularity 
of the collection of random variables {u{t, x)}^t^x)eOT sense of Malliavin. Proposi- 

tion [53] concerns the Malliavin differentiability of order one, while Proposition 15. 71 treats 
higher-order Malliavin derivatives. As already mentioned, the two results rely on hy- 
potheses (a) and (b) of Theorem 15.11 respectively, but not on the lower bound for the 
stochastic convolution. 

Proposition 5.5. Assume that 

(i) no G C(0); 

(ii) / is increasing and belongs to CpQj(M); 

(iii) Q is positivity preserving. 

Let u G C\q>i^q ^6 ^he mild solution to ()2.ip . Then, for any {t,x) € Ot, one has 
u(t,x) € D^'°°. Moreover, the Malliavin derivative Du{t,x) satisfies the following linear 
equation in H: 

Du{t,x) = vo{t,x) + / Gt-six,y){r] - f'{u{s,y)))Du{s,y)dyds (5.2) 
Jo Jo 

where V(){t,x) := (r, z) i— > Gt~T{x,z) 1[q^^{t), and it holds that, for all q >1, 

sup E\\Du{t,x)fH < oo. (5.3) 
{t,x)eOT 

Proof. Since fx is Lipschitz continuous and of class , slight modifications of the "clas- 
sical" results (cf. e.g. [18]) imply that, for any {t,x) € Ot, ux{t,x) belongs to D^'°° and 
satisfies the following linear deterministic integral equation with random coefficients: 

Dux{t,x) = vo{t,x) + / Gt^s{x,y)iv - fx{u\{s,y)))Dux{s,y)dyds. 
Jo Jo 

Recall that, by Proposition 12.51 one has, for any q > 1, 

E\ux{t,x) -u{t,x)\'^ 
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as A — )■ 0, for all {t,x) £ Ot- Therefore, by Lemma HTTl in order to conclude that u{t,x) 
belongs to B^'°° for all (t, x) G Ot, it suffices to show that for all g > 1 and (t, x) S Ot, 
one has 

supE\\Dux{t,x)\\'jj <oo. (5.4) 

A>0 

Let {/I'^'IfcgN be an orthonormal basis of H, and set 

^p'lit, x) := {Dux{t, x), h^)H, {t, x) G Ot. 

Then ip\{t) := ip^{t,-), < t < T, satisfies the deterministic evolution equation with 
random coefficients 



l^^(t) - A^',{t) + F,{t)^1{t) = ^Ht), ^^(0) = 0, 
where Fx{t) := fx{ux{t, •)) - ?? and <^''{t) := Qh''{t). In fact, one has 



(5.5) 



{vo{t,x),h'') 



H 



{Gt-s{x,-),h^{s))L2 ds 



S{t- s)Qh''{s)ds (x). 



Jo 



Gt-six,y)[Qh''{s)]iy) dy ds 



From now we assume, without loss of generality, that = (if not, it is enough to write 
the corresponding equation for t i— )> e~'^^ipx{t), multiplying by e~^* the term $''(t)). 

Fix CO € Q, and let {s,t) i->- U\(t,s), s < t, denote the family of evolution operators 
generated by the time-dependent linear operator A — F\{t). Then we can write 



Jo 



thus also, denoting the integral kernel of Ux{t, s) (that exists by Proposition I4.3r iv) ) . by 
kx{t,s;-,-), 

rt 



ipiit,x) 



kx{t,s;x,y)<^ {s,y)dyds. 



JO 



This yields 



Duxit,x)\\l = Y,\^iit,x)\^ = Y.\ [ I kx{t,s;x,y)i^\s,y)dyds 



Y,\{kx{t,-,x,-)\^^t]{-),h'')L^ 



km 



i%T\^Ll) I 



\kx{t,s;x,-)\\\2 ds. 



Note that Proposition I4.3( ii) implies kx{t,s) < Gt-s pointwise for all < s < t < T. 
Using that Q is positivity preserving, we are left with 



\Dux{t,x)fH < [ \\Gt-s{x 
Jo 



^■jwl^ds = \\vo{t,x)\\'jj. 



(5.6) 
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Let us now show that ||uo(t, is uniformly bounded over t and x. In fact, one has 

\\vo{t,x)\\l= f \\Gt-r{x,-)\\l, dT= f \\B*Gt-r{x.-)\\ldT 

Jo Jo 



{Gt^r{xr):Be^)ldT= [ Y.{f Gt^r{x,y)[Be\y)dyf dr 

•'^ k>l ''^ 

^ f [S{t-T)Be%xfdT, 



fceN 

where {e'^jfcgr^ is an orthonormal basis of L^. The identities 



E|VFA(t,x)P =e| ^ / [S{t-s)Be'']{x)dwk{s)^ = ^ f [S{t - s)Be'']{xf ds, 



yield \\voit,x)\\'jj = E\WA{t,x)\'^ for all G Or- Thanks to Hypothesis O we infer 

that there exists a constant N = N{q), independent of A, such that 

sup E\\Dux{t,x)\\''^ < N. 

We have thus proved that u{t,x) G 1])^'°^ for all {t,x) G Ot- It is therefore lawful to 
apply the Malliavin derivative to the equation satisfied by u, obtaining 

Du{t,x) = vo{t,x) - / Gt-s{x,y)Df{u{s,y)dyds. 
Jo Jo 

Then, appealing to the chain rule proved in Lemma 14. 2^ we obtain that the Malliavin 
derivative Du{t,x) satisfies equation ()5.2p . 

In order to conclude, we only have to show that estimate (|5.3p holds true. The 
argument is essentially the same as above, hence it is only sketched. Still assuming 
r] = without loss of generality, setting F{t,x) := f'{u{t,x)), one has that ip''{t,x) := 
{Du{t,x),h^)H satisfies the linear deterministic evolution equation with random coeffi- 
cients 

- A/(t) + = ^\t), /(O) = 0. 

Let 9.' (ZVL with P(ri') = 1 such that ^ u{t,x,uj) G C([0,r] x 6) for all uj G 9.'. 

Fix UJ G Vl' . Then (t, x) i->- F{t,x) is positive and continuous, hence bounded on the 
compact set [0, T] x O. One can then construct the evolution operator associated to 
A — F, and proceeding exactly as above one arrives at 

sup E,\\Du{t,x)\fjj < oo, 

(t,x)&OT 

so that the proof is complete. □ 

Remark 5.6. Condition (iii) in Proposition 15.51 above does not need to be considered an 
important restriction. Indeed, this condition is satisfied in the spatially homogeneous 
counterpart when dealing with existence and smoothness of the density for stochastic 
heat and wave equations in M"^ (see e.g. |19j). 



Proposition 5.7. Let m G N, ?n > 2. Assume that 
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(i) no G C{0); 

(ii) / is increasing and belongs to C^j(M); 

(iii) Q is positivity preserving. 

Let u G riggNCq be the unique mild solution to ()2.ip . Then, for any {t, x) G Ot, one has 
u{t, x) G D""'*^ and 

sup ||n(t, x) ||d'".9 < oo > 1. 

(t,x)eC'T 

Proof. Let p G N be such that \f{x) + \ f'{x) \ + • • • + < 1 + \x\p for all r G M. 

By Proposition 15. 5|, we have that u(t,x) G 3^'°° for all {t,x) G Ot- Let us consider 
the regularized equation ()2.4p : since /a needs not have bounded derivatives of order 
two and higher, we cannot apply "classical" results (cf. e.g. [T^) to deduce that, for 
any {t,x) G Ot, u\{t,x) belongs to B™'°°. For this reason, we introduce a further 
regularization: let {C/?}/? be a family of mollifiers as in Lemma 14.51 above. Note that 
fxp ~ f'x* all n > 1, hence f\p has bounded derivatives of every order. Let 

u\p be the unique mild solution in RggNCg to the equation 

duxp{t) - Auxpit) dt + fxi3{ux(3it)) dt = 7]uxp{t) dt + B dW{t), ux/siO) = uq. (5.7) 

We split the rest of the proof in three steps: first we show that, for any (t, x) G Ot, one 
has uxpit, x) u{t, x) in L'' as /3 — )■ 0. Then we obtain the uniform bound 

sup E\\D^ux^{t,x)\\^^„<N, (5.8) 

where is a constant independent of A and /3. Finally we pass to the limit as /3 — )> 
and A 0. 

Step 1. We assume again, without loss of generality, that r/ = 0. It is easily seen that it 
holds 

uxfiit) - ux{t) = / S{t - s){fxiuxis)) - fxi3{uxi3{s))) ds, 
Jo 

hence, recalling that S{t) is contracting in L'^{0) and denoting the norm of this space 

by 11-11, 



\\uxi3{t) - Ux{t)\\ < / \\fx{ux{s)) - fxi3{uxi3{s))\\ds. 

Jo 

This yields, by the triangle inequality, 

\\uxpit) - uxit)\\ < / Wfxpiuxpis)) - fxpiuxis))\\ds + Wfxpiuxis)) - fxiuxis))\\ds 







A 

where 



1 /■* , „ 
-T/ \\uxf3{s) - ux{s)\\ds + 



cT 

I/s ■■= I \\fxi3{ux{s)) - fx{ux{s))\\ds. 



By Gronwall's inequality and some obvious manipulations, one arrives at 



Esup \\ux^{t) -uxit)]]" < e«^/^E/l 

t<T 
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Let US show that E/^ ^ as /? ^ 0: since fx is continuous, fxjs converges uniformly to 

fx as (3 ^ 0. Therefore, as ux{s) G C{0) P-a.s., we also have that the integrand in the 
definition of 1^ converges to zero P-a.s. as /? — )■ 0. Taking into account that 

\\fxpMs))-fx{uxis))\\'<xl + \\ux{sW, 

and that E \\ux{s)\\'' ds < oo, we get, by the dominated convergence theorem, that 
E/^ — > 0, hence also 

Esup||nA/3(t) - uxit)\\'l(cJ^ 

t<T ' 

as /3 — 7> 0, for any g > 1. 

Step 2. For n = 1 it is easily seen that ()5.8p holds true, simply by the previous proposi- 
tion, noting that /^^ = /j^ * C/3 ^ 0. For the sake of clarity, let us explicitly show, in the 
case n = 2, how estimate (|5.8p is implied by the corresponding one with n = 1. Then 
the general induction step will be clear. As before, we shall assume, without loss of 
generality, that rj = 0. Since, as already observed before, /a/3 has bounded derivatives of 
every order, we infer that uxj3{t^x) € D^'°°, the iterated Malliavin derivative D'^uxpit^x) 
takes values in H^'^ and satisfies 

D\xp{t,x)+ I I Gt-s{x,y)f':l^{uxp{s,y)){Duxp{s,y)f^dyds 
Jo Jo 

+ 11 Gt-s{x,y)fxp{uxp{s,y)))D'^uxp{s,y)dyds = Q. 
Jo Jo 

Let {/i^jfcgN be an orthonormal basis of H^"^ and set 

^xpit,x) := {D^uxp{t,x),h'')H'S2, G N. 

Then ^\p{t) := '^\p{t, •) satisfies the following linear deterministic evolution equation 
with random coefficients 

j^xpit) - ^A,3{t) + Fxp{t)Ap{t) = ^%{t), = 0, (5.9) 

where 

Fxp{t) ■■= f'x^Mt,-)), ^i^{t) := f'ip{uxp{t,-mDuxp{tr)r\h'^)m2- 
Then we have 

V%{t) = [ Ux^{t,s)^lf,{s)ds, 
Jo 
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hence also, denoting the kernel of Uxis{t, s) by k{t, s; •, •), 

bA/3(*>^)P = f f kxp{t,s;x,y)^lp{s,y)dyds 
Jo Jo 



< 




kxpit, s; X, ?/)|$^^(s, y)\dy ds^ 



10 Jo 

<{ I I Gt-s{x,y)\^'lfs{s,y)\dydsy 




Jo 

t 



<T 



Jo 



Gt-s{x,y)\^'{p{s,y)\^ dy ds, 



where we have used the estimate kxi^ < G in the first inequality, and Cauchy-Schwarz' 
inequality in the third inequality, recalling that S{t) is contracting in L°° (cf. (j3.ip ). 
Summing over k, Tonelli's theorem yields 

\\D^ux^(.t,x)\\l^, = ^\ip'l^{t,x)f <T [ j Gt^s{x,y)Y,\^%{s,y)\Uyds 

h.^-m Jo Jo 

Gt^s{x,y)\\^xi3{s,y)\\H^2 dy ds. 



JO 
where 

^X(s{t,x) := fip{uxp{t,x)){Duxp{t,x)T\ 
Holder's inequality and Tonelli's theorem then imply 

Ep\A/3(t,x)||V I I Gt-s{x,y)n\^x0{s,y)fjj^,dyds, 

Jo Jo 

that is, 

sup ^\D'^Uxfi{t,x)f^^2 <T sup E||^>A/3(t,x)|||^^2. 
(t,x)£OT {t,x)GOT 

Let us show that the right-hand side is finite: by Cauchy-Schwarz' inequality, we have 

Assume, without loss of generality, A < 1, /5 < f. Since /^'^ = f" * C,j3, and, by Lemma 
14.41 there exists o" G N such that \fx{x)\ ^ 1 + \x\^ , by Lemma 14.51 we also have 
I/a/3(^)I <l + kr- Therefore 

E||$A/3(t,x)||V ^ {l + nuxp{t,x)\^'i''f'\nDuxp{t,x)Cf'\ 

where both terms on the right hand side are uniformly bounded over t, x. A, and /3 by 
results already proved; in fact, as for ux itself, the boundedness of the first term on the 
right-hand side above follows from Proposition 6.2.2 in [Ij, since fxp is also monotone. 
The claim is then verified for n = 2. 

The general case is proved by induction in a completely similar way. In particular, 
assume that, given 3 < n < m, one has the uniform bound 

sup ||liA/3(t,x)||mn-l,q < N, 
(t,x)&OT 
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with independent of A and (3. Let {/I'^lfceN be an ortlionormal basis of H®^, and set 

f%{t,x) := {D''ux0{t,x),h'')H<sn. 

Then ip''{t, •) satisfies an equation of the form (|5.9p . where is a sum of finitely many 
terms depending on uxp and on its Malhavin derivatives of order not greater than n — 1, 
whence 

sup E||$A/3(i,2;)||^«„ < A^, 

{t,x)£OT 

with N a constant that does not depend on A nor on (3. Moreover, by an argument 
completely analogous to one already used before, one shows that 

sup E\\D^ux(3{t,x)\\j^^„ < sup E\\^xi3{t,x))\\j^^^ < N, 

where is the same constant of the previous inequality. 

Step 3. Let q > 1. By the previous steps and Lemma |4.1|, passing to the limit as 
/3 — )■ 0, we obtain ux{t,x) € B"^'"^ for all {t,x) € Ot, and also, by lower semicontinuity 
of the norm with respect to weak convergence, 

E\\D''ux{t,x)\\%„ < liminfE||Z)"nA;3(t,x)||V> 

which implies, together with the last inequality, 

sup E\\D''ux{t,x)\\%,,<N. 

{t,x)£OT 

Recalling that, by Proposition 12.51 ux{t,x) — )■ u{t,x) in L'^ as A — )■ for all G Ot, 

appealing again to Lemma 14.11 we arrive at u{t,x) € B™'"? for all {t,x). Since q is 

arbitrary, we conclude that u{t,x) G B™''°° for all {t,x). □ 

5.2. Analysis of the Malliavin matrix. In this subsection, we shall use a standard 
method in order to prove that the inverse of the Malliavin matrix has moments of all 
orders (see e.g. jl91 Theorem 6.2]). 

Proposition 5.8. Assume that the hypotheses of Theorem I5.ir a) are satisfied, as well 
as condition i5.1\} . Then, for any (t, x) €]0, T] x O, one has 

E\\Du{t, x)\\'^'' <oo yq>l. 

Proof. By [17|. Lemma 2.3.1], it suffices to prove that for any q > 2, there exists £oiQ) > 
such that, for all s < Eq, 

F{\\Du{t,x)fH<e)<e'^. (5.10) 

Let {t,x) €]0,T] X O be fixed. Observe that we are assuming the same hypotheses as 
in Proposition 15.51 Hence, u{t,x) belongs to B^'°° and the Malliavin derivative Du{t,x) 
satisfies equation ()5.2p . Then, using this latter equation, we can infer that, for any 6 > 
sufficiently small, 

\\Du{t,x)\\jj= / II Dru{t,x)\\'j^2 dr > / \\Dru{t,x)\\'j^2 dr > -g{x,6) — I{t,x,6), 
Jo Jt-5 2 
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where g{x,6) is as in assumption (|5.1|) and 



I{t,x,5) := / / Gt-~six,y)i'n - f'{u{s,y)))Dt^ru{s,y)dyds 

Jo Jt~T JO 

Hence, using Chebyshev's inequality, we have, for ah e > 0, 



dr. 



F {\\Du{t,x)\\l < e) <F \ lit,x;5) 



> 



g{x,6) 



sA^-ey'mt,x,6w. 

(5.11) 

Let us now find an upper bound for the p-th moment of I{t, x, 6). For this, we start 
by applying Minkowski and Holder's inequalities, the latter with respect to the measure 
on [t — 5,t] X O given by Gt-six,y)dyds, to obtain that 



E\I{t,x,5)\P < E 



t-5 Jo 



< 



t-5 Jo 



Gt-s{x,y)\v - f'{u{s,y))\\\Dt-.u{s,y)\\L2^[o^s];Ll)dyds 

2p-l 



2p 



Gt-s{x,y) dyds 



X / / Gt-s{x,yM{v + \f'Hs,y))\fP\\Dt^.u{s 



|2p 



dyds 



t~s Jo 
t 



t-5 Jo 



IIL2{[0,5];L^)^ 

Gt_,(x,y)E((7? + |/'(n(s,y))|)2f||A-.n(s,y)||f,(p^,j.^^pdyds, 

(5.12) 



where we have also used the estimate ()3.ip . Thus, applying the Cauchy-Schwarz in- 
equality and appealing to the polynomial growth condition on /' (say |/'(2:)| ^ 1 + l-J^ 
for all z € M), the right hand side of (|5.12p can be estimated, up to a positive constant, 
by 

<52P-1 sup (E||A-.n(.,2/)||f^ )V2 

{s,y)e[t-5,t]xO UU,dJ,^Qj 

X jyi^Gt-s{x,y) (l + {E\u{s,y)f'n'^^) dyds. 
At this point, we can appeal to (|5.3p to get 

sup {nDt-Ms,y)\\%,,,..^.y^' < c{T). 

^s,y)e[t-5,t]xO ^ im'^g) 

Taking into account again estimate ()3.ip , and the uniform bound ()3.2p , we can infer that 

E|/(t,x,<5)r <T S'P. 
Plugging this estimate into ()5.1ip . we obtain that 



\Du{t,x)\\l <e) <T 
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On the other hand, (j5.1|) yields g{x, 5) > CxS'^ . Thus, if we choose S = 6{e, x) sufficiently 
small in such a way that CxS'^ = 4e, we get 



F{\\Du{t,x)\\l<e) <T^^<.,T 



where - — 1 > by hypothesis. Therefore, going back to ()5.10p . it suffices to take 
p = and the proof is completed. □ 

Remark 5.9. We should point out that, in fact, in order to prove the existence of density 
(i.e. Theorem 1 5 . 1 I f a) ) . condition ()5.1|) may be slightly weakened. Precisely, one needs to 
prove that ||Z)u(t, > P-a.s. First, by ()5.6|) and the lower semicontinuity of the 
norm with respect to weak convergence, one gets that for any 6 £ (0, 1) and q > 1, 

Then, similarly as above (see also |191 Thm. 5.2]), one proves that 

E\I{t,x,6)\ <5 f I Gs{x,y)g{y,5) dyds, (5.13) 
Jo Jo 

and that, for any n > 1 and 6 € (0, 1), 

i2 ^ In / g{x,s) 



P(||Z)u(t,rr)||^ < -) < - - E|I(t,rr,5)|. 

Taking limit as n — )■ oo and using (j5.13p . we end up with 

F{\\Du{t,x)\\l = 0)<g{x,6r^5 f [ G,ix,y)g{y, 6) dyds. 

Jo Jo 

In conclusion, u{t, x) has a density provided that the following two conditions are satis- 
fied: 

(a) For any x € O, g{x, 6) > for all 6 > 0, 

(b) For any x € O, it holds 



I™ ~r~xV / / Gsix,y)giy,6) dyds = 0. 
s^o g{x,d) Jq Jo 



5.3. Proof of Theorem 15.11 It is now just a matter of putting pieces together. In 
particular, in view of Bouleau-Hirsch criterion (see e.g. \n\ Thm. 2.1.3]), (a) follows by 
Propositions 15.51 and 15.81 Similarly, (b) follows by Propositions 15.71 and 15.81 applying a 
general criterion of the Malliavin calculus (see e.g. [n\ Prop. 2.1.5] or \13\ Thm. 4.1]). 
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